In order to realize the main physics goal at the future hadronic colliders 
LHC and SSC, the precise estimate of inclusive production of neutral clusters 
is needed to pin down signals due to Higgs particle or new physics Q. For this 
purpose leading order (hereafter denoted as LO ) perturbative QCD predictions 
- based on evaluations of partonic cross sections at tree level and evolution 
of structure and fragmentation functions at one loop level- are not accurate 
enough. A consistent calculation at next to leading order (hereafter denoted 
as NLO) needs two loop evolved structure and fragmentation functions and a 
NLO evaluation of parton-parton subprocesses which has been performed by 
;h . Aversa et al. [|| a few years ago and is based on the 0(af ) matrix elements of 

Ellis and Sexton [|. 

, Up to now such an analysis has been performed for for heavy quarks and 

' pion fragmentation functions ||5|. Following the latter approach we present a 

complete NLO evaluation of eta inclusive production at hadronic colliders in 
order to estimate, as precisely as possible, the rj rates at LHC and SSC. 

We first extract ry fragmentation functions from e~^e~ and pp collisions at the 
' scale =30 GeV using one of the three methods developed in our previous 

■ work|5|, where we have shown that those different approaches give consistent 

f-f^ I theoretical predictions for vr*^ production at LHC/SSC. Therefore we are quite 

0^ ' confident on the accuracy of the procedure followed below for rj. 

r"| ' Then we will confront our predictions with the existing experimental data in 

P^, e'*"e~ and pp annihilation in order to estimate the theoretical uncertainty. Cross- 

Q^l sections for r] production in hadron collisions are usually not directly measured, 

^ ■ but they are deduced from vr*^ cross-sections, assuming an experimental rj/n'^ 

ratio which has been measured at ISR@ and is taken to be constant of order 
0.5. 

^ ■ We find that the theoretical rj/ir'^ ratio increases with pt at present energies, 

and therefore we study in detail this ratio in the LHC/SSC energy range, when 
we give our predictions for the future colliders. 

We briefly report, for reader's convenience, the main formulae for one hadron 
inclusive production at next-to-leading-order|5|, via the generic reaction A + 
B — > H+X, where A and B stand for hadrons and/or leptons. The cross section 
is given by the convolution of the partonic cross section and the fragmentation 
functions Dl^{z,Mj): 

(1) 

where zu is the reduced energy of the hadron H, zh = lEuj^fS^ 9 is the 
scattered angle of the parton 1, and the inclusive production of the parton 1 in 
the reaction A + B ^ I has the following perturbative development: 

i^.%|^(^,^,«.(^^),M|,..0=aVB^,(^,^)+^4^B-..(^,^,M|)+. 

(2) 
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Furthermore Dj^ [z, Mj) represents the number of hadrons H inside the parton 
1 carrying the fraction of momentum z from H, evolved at the scale Mj. These 
fragmentation functions satisfy Altarelli-Parisi type evolution equations 



For inclusive production in e^e annihilation, the partonic cross-sections at 
next-to-leading order are given by: 
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(3) 



(4) 



where o"o is the usual point like cross-section 



a is the QED coupling constant and is the invariant mass of the e^e~ pair. 
The functions KJ, K^, KJ and Kg have been extracted from the reference Q 
(see also HI) and are given elsewhere||5|. 



For hadronic collisions the partonic cross-sections are: 
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where the partonic variables are s = X1X2S and 

1 - V 



VW 

Xl = , X2 

VW 

and the hadronic ones are defined by: 



1-v' 
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2^ ^' 2 -2/(1 -cos 6*) 
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We now consider the inclusive production in e'^e annihilation at MfQ = 
V5 = 30 GeV, as simulated by the Monte Carlo generator HERWIGg. As 
well known, this event generator includes the QCD parton shower to leading 
and next to leading accuracy - in particular the kinematical corrections due to 
the phase space boundaries are summed up to all orders - as well as the hadro- 
nisation of the color singlet clusters into the physical particles. Furthermore 
HERWIG has been shown |1C] to describe with good accuracy the observed 



features of PETRA and LEP data. Then, similarly to what have been done for 
7r° inclusive production, we will use the rj distribution generated by each quark 
flavor which originates from the photonic vertex, as a realistic description of the 
quark fragmentation into r/. Owing to the symmetry of quarks and antiquarks 
fragmenting into r] we extract the quark fragmentation functions from: 

'^^^^g-l^{zH, M%) ~ 6ao E elD^^izH. M%), (6) 

where o"o is the pointlike cross section defined above. 

The reaction e^e~ — > r/ + X has been therefore decomposed into each 
contribution e"^e~ uu, dd, ss, cc and bb. The generated distributions are 
parametrized as 

and analyzed using the minimization procedure MINUIT. The coefficients Ni 
are constrained by the normalization condition: 

^'^dzDi{z,M%) = {n^)„ (8) 

where the average values (ra^)^ are given by HERWIG for each quark flavor, in 
agreement with the total observed multiplicity (n-j^). The parameters Ni, Ui and 
f3i are extracted from the r] inclusive distribution generated, for each flavor, in 
the z range .025 < Zfj < -95 and shown in Table I. As can be inferred from this 
table the statistical error on the parameters is less than 5%. 



D'l{z,M%) = N,z''^{l-zt (7) 



So far we have not included the contribution from gluon. Indeed from the 
analysis of the three jet events it would be possible, in principle, to extract 
from HERWIG the appropriate information. The corresponding accuracy is 
however unsatisfactory, due to the limited sensitivity to hard gluon effects in 
e^e~ annihilation. Then, as for the tt*^ analysis, we have followed a different 
approach to extract from HERWIG the gluon fragmentation function. We have 
analyzed the subprocess gg ^ gg ^ r] + X from pp annihilation at MfQ = ^/s ~ 
30 GeV, in analogy to the quark case. In order to eliminate the background 
from the fragmentation of the spectator partons we have constrained the etas 
to lye within a cone of semi aperture 5 = .35 — .40 rad around the direction of 
the parent gluons emitted at 90 deg. The value of 5 is found by an appropriate 
angular study of the generated distribution. With a parametrization of the 
form (0) we find the values of the parameters Ng, ag and Pg given in Table II, 
which define Sets I and II of fragmentation functions. 



3 



After inclusion of the gluon fragmentation function and using NLO evolution 
together with NLO terms in the rj inclusive cross section (eqs. 3, 4) we show in 
Figure 1 our results at = 35 GeV, compared with JADEQ and CELLO 



12(1 data. The agreement is satisfactory as can be inferred from the figure. The 



difference between Sets I and II is negligeable . 

After evolution to \/5 = 91.2 GeV, we also obtain good agreement with L3 
1 13] LEP data as shown in Fig 2. 



We consider now our predictions for inclusive production in hadronic col- 
liders. We first compare with data from CERN ISR §, for ^/^ = 52.7 GeV 
and -k/S = 62.4 GeV, as shown in Figures 3 and 4 for /i = M = Mf = Pt and 
fi = M = Mf = Pt/2 using the quark fragmentation functions from Table I and 
the two gluon solutions from Table II, for 6 = 0.35 (Set I) and 5 = 0.40 (Set 
II). In doing so, in absence of direct data on t] production, we have inferred 
the cross section from vr*^ data assuming the experimental]^ ti/tt^ ratio R of 
0.58 lb 0.05 and 0.55 it 0.06 respectively, indipendent from pt- The agreement is 
satisfactory within the theoretical and experimental uncertainties. 

Let us focus now on the UA2 data at the SppS collider |l^]. We will use 
two sets of quite precise vr'^ data, for Pt < 15 GeV and pseudorapidity y ~ 
and for 15 < < 45 GeV and y ~ 1.4, and an experimental ratio 'q/ir^ of 0.5 
as obtained from ISR dataQ. The comparison with the theoretical predictions 
is shown in Figures 5 and 6 for fj, = M = Mf = Pt/2, Pt and for the two gluon 
sets of fragmentation functions. The agreement is quite good at low pt, and 
slightly favours set I. 

On the other hand we note that at higher pt the comparison with data 
suggests a larger value for the ratio rj/ir^, and therefore apt dependence for this 
ratio. Indeed from the result of our previous study on inclusive vr'^ production, 
we show in Figure 7 the predicted pt dependence of P = 77/7r at =630 GeV, 
which indeed rises with pt . 

Finally we proceed to the predictions for LHC and SSC at ^/S = 16 TeV 
and = 40 TeV respectively. The cross sections are calculated at LO (Born) 
and NLO and using HMRS Set of structure functions [|^] and are displayed in 
Fig. 8.|3 

To estimate the theoretical uncertainty we study in Fig. 9 the ratio of the 
two predictions from the two different choices of gluon fragmentation functions, 
evolved to NLO accuracy, at \^=16 TeV. 

We then show the theoretical ratio rj/ir^ as evaluated from the vr'^ results 
of reference |^. As one can infer from Fig. 10 the ratio increases with respect 
to ISR energies, and shows a dependence on pt similar to what found at SppS 
energies. 

The uncertainty due to factorisation scheme, especially coming from frag- 
mentation functions is expected to be tiny because the evaluation done for one 

^The discontinuities in the curves are simply due to CPU time limitation on the number 
of the data points. 
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jet inclusive cross section has shownQ that at coUider energies its magnitude is 
of the order of 5% and we can reasonably expect the same order of magnitude 
for one hadron inclusive cross section. Finally, the theoretical uncertainty from 
the structure functions is much smaller than that coming from fragmentation 
function. 

To conclude, we have performed a complete next to leading order analysis of 
inclusive r] production in e~^e~ and hadronic collisions. We have found that the 
theoretical ratio r]/TT^ depends on pt and we have presented results for the future 
colliders, where the absolute rates can be predicted within a factor of two. This 
will certainly be of help for neutral background rejection at supercolliders. 
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Table captions 

• Table I: parameters of the quark fragmentation functions as obtained from 
HERWIG in e+e" annihilation at Mq = 30 GeV. 

• Table II: parameters of the gluon fragmentation functions as obtained 
from HERWIG at Mq = 30 GeV, with two hypotheses on the angle 6 (see 
text). 

Figure Captions 

• Fig. 1: NLO inclusive rj production in e'^e" annihilation with the quark 
and ghion fragmentation functions evolved at s/S = 35 GeV, compared 
with data. 

• Fig. 2: NLO inclusive rj production in e~^e~ annihilation with the quark 
and gluon fragmentation functions evolved at \/S = 91.2 GeV, compared 
with data. 

• Figs. 3: NLO inclusive r] production in pp collisions at ISR energies for 
H = M = Mf = Pt,Pt/2 for a) 6 = 0.35 (Setl) and 6 = 0.40 (Set II), see 
text. 

• Figs. 4: NLO inclusive t] production in pp collisions at ISR energies for 
H = M = Mf = Pt, Pt/2 for a) Set I and b) Set II. 

• Figs. 5: NLO inclusive rj production in pp collisions at SppS energies for 
^ = M = Mf = Pt, Pt/2 for Sets I and II, at ^/S = 540 GeV and y = 0. 

• Figs. 6: same as figs 5 at VS = 630 GeV and y = 1.4. 

• Fig. 7: theoretical prediction of the ratio rj/ir'^ as function oi pr at SppS 
energy. 

• Fig. 8: NLO inclusive rj production in hadronic collisions at LHC and 

SSC energies at y = 0. 

• Fig. 9: ratio of inclusive rj cross sections for the two sets of gluon frag- 
mentation functions. 

• Fig. 10: theoretical prediction for the ratio tj/tt^ at LHC energy. 
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Process 


a 


/3 




< > 


e+e^ 


uu 


-0.91 ±0.02 


2.09 ±0.07 


0.24 


0.35 


e~^e~ — > 


dd 


-0.88 ±0.02 


2.14 ±0.08 


0.26 


0.36 


e"'"e~ 


ss 


-0.72 ±0.02 


2.73 ± 0.08 


0.37 


0.42 


e+e~ 


cc 


0.14 ±0.03 


7.10 ±0.14 


8.73 


0.57 


e+e~ — > 


bb 


-0.20 ±0.05 


11.24 ±0.31 


9.92 


0.69 



Table I 
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a 






< > 


I 


0.35 rad 


-0.18 ±0.06 


4.58 ± 0.25 


2.52 


0.51 


II 


0.4 rad 


-0.43 ± 0.06 


3.47 ± 0.26 


1.48 


0.62 



Table II 
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